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Abstract 

Water  transport  through  the  membrane  of  a  PEMFC  was  investigated  by  measurement  of  the  net  drag  under  various  feed  gas  humidity.  Measured 
data  were  compared  with  computed  results  obtained  using  a  two-dimensional  cell  model.  Considering  the  change  in  the  gas  content  related  to  the 
flow  configuration,  the  humidity  of  the  supply  gas,  reaction  rates,  and  the  mass  balance  of  each  gas  species  were  derived  at  five  sections  along  the 
flow  channels.  By  solving  these  mass  balance  equations,  the  water  transport  rates  and  current  density  distribution  were  obtained  along  the  flow 
channels  for  various  feed  gas  humidity.  The  results  for  net  drag  computed  from  the  model  show  a  similar  tendency,  but  are  slightly  higher  than  the 
measured  values.  This  suggests  that  there  is  a  certain  resistance  related  to  water  transport  at  the  cathode  membrane  interface  in  association  with 
water  production.  The  cause  of  this  water  transport  resistance  is  discussed. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  has 
attracted  a  great  deal  of  attention  in  the  last  decade  as  a  promis¬ 
ing  candidate  for  a  high-efficiency  low-emission  power  source 
in  both  mobile  and  stationary  applications.  For  mobile  use,  fuel 
cells  are  usually  operated  under  a  relatively  high  current  den¬ 
sity  in  order  to  achieve  higher  power  density.  For  stationary  use, 
fuels  cell  are  usually  operated  under  moderate  current  density,  in 
order  to  achieve  higher  efficiency.  However,  there  are  still  many 
problematic  issues  that  must  be  overcome.  One  major  problem 
is  the  water  management  of  fuel  cells.  The  operation  of  a  state- 
of-the-art  PEMFC  requires  careful  water  management.  At  low 
humidity,  the  proton  exchange  membrane  and  electrode  assem¬ 
bly  (MEA)  lose  water,  which  leads  to  a  rapid  increase  in  ohmic 
resistance.  Conversely,  if  too  much  liquid  water  is  present  in 
the  cell,  the  pores  in  the  electrodes  will  be  filled  with  water 
and  the  passage  of  reactant  gases  obstructed.  In  other  words, 
the  operating  conditions  and  MEA  components  have  to  be  well 
matched  in  order  to  avoid  membrane  dehydration  and  cathode 
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flooding.  In  addition  to  the  experimental  approach,  models  for 
water  management  are  useful  to  understand  the  processes  that 
govern  water  transport.  A  clear  understanding  of  water  trans¬ 
port  processes  will  assist  the  optimization  of  fuel  cell  operating 
conditions  and  relevant  electrode  structure. 

Several  experimental  results  have  been  reported  in  the  lit¬ 
erature  for  the  water  transport  in  proton  exchange  membranes 
(PEM),  such  as  the  net  drag  coefficients  measured  by  Janssen 
et  al.  [1].  Although  their  investigation  was  mainly  related  to 
the  effect  of  humidity,  the  humidity  of  the  inlet  gas  was  limited 
to  a  dry  or  wet  condition.  Ren  et  al.  reported  measurements 
for  the  electro-osmotic  drag  of  water  of  a  PEM  in  a  direct 
methanol  fuel  cell  (DMFC)  [2].  The  experiments  were  per¬ 
formed  under  the  condition  that  the  water  flux  across  the  mem¬ 
brane  was  exclusively  driven  by  the  electro-osmotic  drag.  Dong 
et  al.  reported  current  distribution,  species  distribution,  and  HFR 
(high  frequency  resistance)  data,  using  a  single  serpentine  flow 
channel  in  a  co-flow  arrangement,  under  mainly  low-humidity 
conditions  [3].  Experimental  data  suggest  that  humidification 
of  the  anode,  rather  than  the  cathode,  was  demonstrated  as 
a  most  critical  factor  for  achieving  high  performance  at  the 
inlet  region. 

It  is  useful  that  these  experimental  results  are  compared 
with  computed  results  from  a  cell  model.  In  addition  to  exper- 
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Nomenclature 

aj  activity  of  water  in  stream  j  (anode,  cathode) 

C  concentration  of  water  in  the  membrane 

(mol  cm-3) 

Dw  diffusion  coefficient  of  water  (cm2  s-1) 
e  ratio  of  actual  to  apparent  surface  area 

fj  fraction  j  (solid  s,  open  o)  of  porous  surface 

F  Faraday  constant  (96,485  C  mol- 1 ) 

io  exchange  current  density  (A  cm-2) 

I  local  current  density  (A  cm-2) 

m  molar  flow  rate  (along  flow  direction)( 

mol  s-1  cm-2) 

rid  electro-osmotic  drag  coefficient 

Nw  molar  flux  of  water  (perpendicular  to  the  MEA) 

(mols-1  cm-2) 
p  pressure  (kPa) 

r  pore  radius 

R  gas  constant  (8.314  J  mol-1  K-1) 

S  surface  area  of  the  catalyst  layer  (cm2) 

td  dew  point  of  feed  gas  (°C) 

T  temperature  (K) 

Vo  open  circuit  voltage  (V) 

v  coordinate  along  the  gas  channel 

y  coordinate  perpendicular  to  the  MEA 

Greek  symbols 

a  net  drag  of  water  per  proton 

y  surface  tension 

8  thickness 

£  porosity  of  the  electrode 

0  over-potential  for  the  ORR,  V 

f  stoichiometry 

9  contact  angle 

p  density  of  the  catalyst  layer 

cr  conductivity  of  the  membrane  (£2~l  cm-1) 

r  tortuosity  factor 


Subscripts 

a 

anode 

c 

cathode 

in 

inlet 

m 

membrane 

out 

exit 

sat 

saturation 

V 

vapor 

Nguyen  et  al.  developed  a  combined  heat  and  mass  transfer 
model  for  a  co-flow  configuration  under  the  condition  that  all  of 
the  water  in  the  electrodes  is  in  the  vapor  phase  [7].  Janssen  pre¬ 
sented  a  two-phase,  two-dimensional,  along-the-channel  model 
[8].  The  model  explains  the  water  transport  in  electrodes  under 
saturated  and  non-saturated  gas  conditions,  by  applying  the 
concentrated  solution  theory  under  the  assumption  of  uniform 
current  density  at  all  locations.  During  PEMFC  operation  with 
unsaturated  reactant  gases,  especially  at  low  stoichiometric  flow 
rates,  liquid  water  is  likely  to  appear  in  the  cathode.  Pasaogullari 
and  Wang  modeled  the  two-phase  flow  and  transport  in  the  air 
cathode  of  a  PEMFC,  and  concluded  that  capillary  action  is  dom¬ 
inant  inside  the  two-phase  zone  [9] .  Recently,  three-dimensional 
models  have  been  developed  to  investigate  the  distribution  of 
current  density  and  temperature  in  the  MEA  and  gas  diffusion 
layer  (GDL),  and  the  distribution  of  gas  velocity  and  gaseous 
species  in  the  cross  section  of  the  GDL  and  flow  channel  under 
higher  power  density  [10,1 1]. 

There  have  been  few  reports  concerning  the  comparison  of 
water  net  drag  experimental  results  and  computed  simulation 
results.  The  purpose  of  the  present  work  is  to  compare  experi¬ 
mental  net  drag  data,  obtained  under  various  reactant  gas  humid¬ 
ity  and  different  operation  modes,  with  the  results  computed 
using  a  relatively  simple  two-dimensional  cell  model. 

2.  Measurements  of  net  drag 

2.7.  Experimental  apparatus  for  net  drag  measurement 

The  change  of  water  content  in  both  the  anode  and  the  cath¬ 
ode  flow  channels  was  measured  under  various  cell-operating 
conditions.  A  25  cm2  MEA  was  constructed  using  a  Nation  1 12 
membrane  with  a  platinum  loading  of  0.5  mg  cm-2  on  both 
the  anode  and  cathode.  The  MEA  was  mounted  between  two 
graphite  current  collector  plates,  each  of  which  contained  a  ser¬ 
pentine  channel  machined  to  a  width  of  1  mm  and  a  depth  of 
1 .0  mm.  The  MEA  was  installed  between  two  gaskets  in  order 
to  prevent  gas  leakage.  A  counter  flow  configuration  was  used. 
Table  1  gives  the  specifications  for  the  MEA. 

The  cell  was  operated  at  constant  temperature,  which  was 
controlled  by  an  electric  heater.  The  feed  gas  humidity  was  con¬ 
trolled  by  bubbling  of  the  gas  through  water  that  was  maintained 
at  a  set-point  temperature.  At  the  exits  of  the  gas  flow  channels, 
the  water  vapor  was  condensed  in  cold  traps  and  the  amount  of 
water  was  measured  at  both  the  anode  and  the  cathode.  Before 
loading  the  platinum  catalyst,  a  thin,  hydrophobic,  microporous 
layer  was  applied  to  the  interface  of  the  catalyst  layer/GDL. 


imental  investigation,  much  research  has  been  performed  for 
PEMFC  modeling,  which  range  from  one-dimensional  to  three- 
dimensional  models  [4-11].  A  one-dimensional  model  neglects 
the  variation  in  the  water  content  of  the  membrane,  and  the  gas 
concentration  along  the  flow  channel.  In  order  to  investigate  the 
distribution  of  gaseous  species,  the  water  content  of  the  mem¬ 
brane,  the  current  density  and  the  net  water  transport  along  the 
flow  channel,  two-dimensional  models  have  been  developed. 


Table  1 

Parameters  used  for  the  test  cell 


Parameter 

Symbol 

Value 

Membrane  equivalent  weight 

E 

1200 

Membrane  thickness 

50  x  10-4  cm 

Membrane  dry  density 

P 

1.84 

Catalyst  layer  thickness 

^cat 

10  x 10-4  cm 

GDL  thickness 

GDL 

300  x  10-4  cm 

Tortousity  factor 

T 

5 
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Cathode  Humidified  Temperature  (°C) 


Cathode  Humidified  Temperature  (°C) 


Fig.  1.  Net  water  drag  a,  as  a  function  of  cathode  humidified  temperature  (FC 
mode,  /cell  =  80  °C,  /da  =  60  °C,  £f  =  1.43,  £a  =  2);  (□)  calculation;  (•,  O)  exper¬ 
imental  results;  (•)  derived  from  the  water  balance  at  the  cathode;  (O)  derived 
from  the  water  balance  at  the  anode. 


Fig.  2.  Net  water  drag  a,  as  a  function  of  cathode  humidified  temperature  (H2 
mode,  /Ceii  =  80  °C,  /da  =  ^dc>  £f=1.43);  (□)  calculation;  (•,  O)  experimental 
results;  (•)  derived  from  the  water  balance  at  the  cathode;  (O)  derived  from 
the  water  balance  at  the  anode. 


In  order  to  determine  the  water  transport  mechanism,  two 
types  of  experiments  (FC  and  H2  mode)  were  performed.  In  H2 
mode  operation,  humidified  hydrogen  and  nitrogen  were  sup¬ 
plied  as  anode  and  cathode  feed  gases,  using  the  same  bubbling 
water  temperature,  and  a  direct  voltage  was  applied  between  the 
two  electrodes.  In  this  mode,  no  generation  of  water  occurred  at 
the  cathode.  Therefore,  flooding  in  the  GDL  was  not  expected  to 
have  occurred  in  this  operation  mode.  In  the  FC  mode,  humidi¬ 
fied  hydrogen  and  air  were  supplied  as  in  conventional  fuel  cell 
operation.  Measurements  were  conducted  for  12  hours  under 
each  of  the  specified  conditions. 

2.2.  Experimental  results 


Net  water  drag  was  derived  from  the  measured  water  content 
for  both  FC  and  H2  mode  cell  operating  conditions  at  various 
feed  gas  humidity.  The  molar  flux  of  water  is  given  by  the  fol¬ 
lowing  equations  for  the  cathode  and  anode,  respectively: 


I  I 

Nw  c  —  —  a  — 

w,c  2F  F 


(la) 


al 

~F 


(lb) 


where  the  first  term  in  Eq.  (la)  is  the  flux  of  water  generated  in 
the  cell  reaction  and  the  second  term  is  that  of  net  drag  transport. 
Because  the  total  molar  flux  of  water  are  equivalent  to  the  dif¬ 
ference  of  water  contents  at  the  exit  and  inlet,  the  net  drags  are 
derived  for  the  cathode  and  anode  water  contents,  respectively: 


F  1 


(mH20,c,ou,  «H2O,0,in)  2 

(lc) 

(^H20,a,in  ~  ^H20,a,out) 

(Id) 

At  a  current  density  of  0.3  A  cm-2,  the  stoichiometry  at  the 
anode  and  cathode  are  1.43  and  2,  respectively.  The  values  of 
net  drag  are  derived  from  the  collected  water  weight  at  both 
the  cathode  and  anode  exits.  The  balance  of  the  water  collected 
at  the  anode  and  cathode  was  between  98.3  and  108.5%  for 
various  feed  gas  humidity.  Figs.  1  and  2  show  the  net  drag 
measured  at  a  cell  operating  temperature  of  80  °C  during 


FC  and  H2  mode  operation,  respectively.  The  range  of  each 
measured  value  indicates  the  difference  between  the  anode  and 
cathode  data,  derived  from  the  change  in  the  amount  of  water 
at  the  inlet  and  the  exit.  Computed  results,  based  on  the  cell 
model,  are  also  provided  in  these  figures.  The  cell  model  is 
described  in  the  following  section. 

For  FC  mode  operation,  the  anode  gas  humidity  was  main¬ 
tained  at  tda  =  60°C,  while  the  cathode  gas  humidity  ranges 
between  60  and  75  °C.  The  humidity  is  expressed  using  the  dew 
point  of  the  inlet  cathode  (anode)  gas  as  t&c  (/da)-  The  values  of 
net  drag  decreased  as  the  cathode  gas  humidity  increased.  For  a 
humidity  higher  than  tdc  =  75  °C,  the  net  drag  has  negative  val¬ 
ues,  which  means  water  transfers  from  the  cathode  to  the  anode. 

During  H2  mode  operation,  the  humidity  of  the  cathode  gas 
was  changed  under  the  condition  /da  =  Gc-  Contrary  to  FC  mode 
operation,  the  values  of  net  drag  increased  as  the  humidity  of 
the  feed  gas  increased. 

3.  Discussion 

3.1.  Cell  model 

In  order  to  evaluate  the  experimental  results,  a  relatively  sim¬ 
ple  two-dimensional  cell  model  was  developed.  The  schematic 
model  is  shown  in  Fig.  3.  The  model  consists  of  a  membrane 

anode  flow  cathode  flow 

channel  membrane  channel 


Fig.  3.  Schematic  model  of  PEMFC.  x  is  a  coordinate  along  the  gas  channel  and 
_y  is  a  coordinate  perpendicular  to  the  membrane. 
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sandwiched  between  two  gas  diffusion  electrodes  (catalyst  lay¬ 
ers  and  gas  diffusion  layers)  and  two  flow  channels.  The  model 
accounts  for  the  mass  transport  of  water  through  the  membrane 
(y-direction),  the  water  content  in  the  membrane  and  the  current 
density  distribution  along  the  flow  channel  (v-direction).  The 
baseline  reactant  gases  of  the  anode  and  cathode  are  consid¬ 
ered  to  be  humidified  hydrogen  and  humidified  air  or  nitrogen, 
respectively. 


3.1.1.  Assumptions 

In  order  to  keep  the  water  transport  model  simple,  a  number 

of  assumptions  are  made: 

1.  The  temperature  is  assumed  to  be  uniform  throughout  the 
ME  A  and  the  flow  channels.  This  assumption  is  rationalized 
for  a  small  single  cell  maintained  at  constant  temperature  by 
an  electrical  heater. 

2.  An  ideal  gas  mixture  is  assumed. 

3.  The  water  content  in  the  membrane,  at  the  interfaces  of  both 
electrodes  is  assumed  to  be  in  equilibrium  with  water  vapor 
in  the  gases. 

4.  The  pores  of  the  catalyst  layer  and  the  GDL  for  gas  diffusion 
are  separated  from  the  liquid  flow  passages.  The  existence 
of  liquid  water  only  modifies  the  effective  porosity  for  gas 
diffusion. 

5.  The  gases  and  liquid  water  are  treated  as  being  perfectly 
mixed  and  are  assumed  to  flow  with  the  same  velocity. 


3.1.2.  Water  transport 

Water  is  transported  in  the  membrane  by  electro-osmotic  drag 
caused  by  proton  transport  and  back  diffusion  due  to  the  con¬ 
centration  gradient  of  water. 

_  /(*)  n  dC(x)  _  I(x) 

N w,a  —  _  D  w  —  a  (2) 

F  oy  F 

where  a  is  the  net  drag  coefficient. 

The  electro-osmotic  drag  coefficient  (nf)  and  the  diffusion 
coefficient  ( Dw )  are  given  by  Springer  et  al.  [4],  using  the  activity 
of  water,  aj,  where  j  is  applicable  to  both  the  anode  (a)  and  the 
cathode  (c)  as  follows: 


rid  =  0.0049  +  2.024 aj  —  4.53  aj  +  4.09  Oj, 
=  1.59  +  0.159 (aj  -  1), 


aj  <  1 

Clj  >  1 


aj  —  Cj 


P 

i+at 


_ H 

Dw  =  ndj  x  5.51  x  10 


x  exp 


(2416  x 


(4) 

(5) 


/,  in  the  cathode  (direction  of  flux  is  defined  from  anode 
to  cathode): 


Nq2  (x) 


m 

4  F 


^W,cC*0 


IF 


(1  +  2a) 


For  a  PEM  electrode,  Sp-21.2  x  104 cm-1  and 
e  =  0.27-0.35  are  reported  [12].  In  the  catalyst  layer,  micro 
carbon  support  particles  aggregate  to  form  the  agglomerated 
structure,  the  peak  diameter  of  which  is  approximately  0.07  [xm, 
derived  from  an  experimental  pore  distribution  curve  [12,13]. 
Since  the  mean  free  path  of  H2O  at  101.3  kPa  is  approximately 
0.1  pan,  it  is  reasonable  to  consider  Knudsen  diffusion.  On 
the  other  hand,  for  the  GDL,  it  is  sufficient  to  consider  only 
molecular  diffusion,  because  the  mean  diameter  of  the  GDL 
pores  is  approximately  30  pan.  In  an  operating  cell,  the  porosity 
of  the  cathode  catalyst  layer  is  expected  to  be  less  than  that 
determined  by  BET  (Brunner-Emmet-Teller)  measurements, 
because  a  number  of  pores  are  filled  with  liquid  water  produced 
during  the  cell  reaction.  Assuming  that  £  =  0.1,  the  estimated 
^eff,H2o,o2’  A?ff,o2,N2  and  De ff,H2o,H2  are  0.0043,  0.0033  and 
0.0134  cm2  s-1,  respectively,  which  is  less  than  1/60  of  the 
molecular  diffusion  coefficients. 

The  effective  diffusion  distance  in  the  direction  of  thickness 
is  expressed  as  <$gdl  x  t,  using  the  tortuosity  factor  (r).  Var¬ 
ious  values  of  the  tortuosity  factor,  between  2.5  and  7,  have 
been  reported  [1,14].  However,  for  the  cell  model,  these  values 
should  be  modified  considering  the  cell  configuration.  Usually, 
in  a  PEMFC,  reactant  gases  are  supplied  through  gas  channels 
that  are  constructed  using  machined  carbon  current  collector 
plates  adjacent  to  the  GDL.  Therefore,  at  the  GDL/channel  inter¬ 
face,  there  exists  a  channel  section  that  contacts  the  reactant 
gas,  and  a  collector  section  to  which  the  reactant  gas  is  not 
directly  supplied.  The  reactant  gases  supplied  from  the  channel 
section  are  transferred  to  the  membrane/electrode  interfaces,  dif¬ 
fusing  toward  the  rear  of  the  collector  section.  The  gas  diffusion 
effect  in  the  planar  direction  can  be  estimated  by  comparing 
the  computed  oxygen  concentration  at  the  GDL/catalyst  layer 
using  the  two-dimensional  (y-z-direction)  diffusion  equation, 
and  that  computed  by  the  one-dimensional  (y-direction)  dif¬ 
fusion  equation.  The  difference  between  these  concentrations 
is  equivalent  to  the  difference  between  the  effective  diffusion 
distances.  The  increase  in  the  effective  diffusion  distance  is  esti¬ 
mated  to  be  approximately  15%.  Based  on  an  average  of  the 
reported  tortuosity  factors,  and  the  increase  in  effective  diffusion 
distance,  a  tortuosity  factor  of  r  =  5  is  applied  to  the  following 
computation. 


3.1.3.  Gas  diffusion 

Gas  diffusion  in  the  y-direction  is  expressed  by  the 
Stefan-Maxwell  equation.  In  the  cathode,  the  diffusion  equa¬ 
tion  for  the  ternary  gas  mixture  is  simplified  under  the  condition 
of  flux  Vn2  =  0  (N2  is  stagnant).  In  the  anode,  the  diffusion 
equation  for  a  binary  gas  mixture  is  applied.  The  following 
relationships  exist  between  Nq2  and  Ajpo  and  current  density 


3.1.4.  Distribution  of  current  density  and  water  transport 
The  relationships  shown  in  Eqs.  (8)-(10)  exist  between  the 
cell  output  voltage  V,  the  current  density  /,  the  cell  over  potential 
r/,  and  the  membrane  ion  conductivity  a.  The  distribution  of  the 
current  density  along  the  flow  channel  is  calculated  using  Eqs. 
(8)-(10)  under  the  condition  of  a  uniform  cell  voltage  along  the 
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flow  channel  and  a  given  total  current: 


V  =  Vo- o- 


ftn7(v) 

a 


a  =  (0.00514  x  8.8 nd  -  0.00326) 
X«XP  (l267  x  (jij  -  1)) 


2RT  (  I{x) 
ii  =  - In  - 

F  \iopo2 


(10) 


In  the  cell  model,  the  cell  is  divided  into  five  sections  in 
the  v-direction.  In  each  section,  the  mass  balance  of  each  gas 
species  is  expressed  as  a  set  of  non-linear  equations.  These  sets 
of  equations  are  solved  by  the  Newton  method  to  obtain  the 
distribution  of  current  density,  net  water  drag  and  the  contents 
of  gas  species  along  the  flow  channel. 


3.2.  Comparison  of  the  model  computation  and 
experimental  measurements 

3.2.1.  Water  drag  dependence  on  humidification 
temperature 

Comparison  of  the  measured  net  drag  with  the  computed 
results  derived  from  the  previously  described  cell  model  are 
shown  in  Figs.  1  and  2.  Although  the  calculated  results  are 
slightly  higher  than  the  measured  values,  they  do  demonstrate 
the  tendency  well,  especially  for  the  H2  mode.  The  discrepancy 
in  the  computed  results  from  the  FC  mode  is  larger  than  that 
from  the  H2  mode.  This  suggests  that  there  is  a  resistance  of 
some  kind  for  the  transport  of  generated  water  to  the  cathode 
channel  at  the  cathode  electrode/membrane  interface. 

The  pressure  of  the  water  in  the  catalyst  layer  is  higher  than 
the  saturation  pressure,  due  to  the  existence  of  the  hydrophobic 
microlayer  at  the  cathode/membrane  interface  (capillary  effect). 
This  means  that  the  activity  of  water  at  the  cathode/membrane 
interface  is  increased  by  the  existence  of  this  hydrophobic  micro¬ 
layer.  The  pressure  of  the  water  is  expressed  in  Eq.  (11),  using 
surface  tension  (y),  and  apparent  contact  angle  (ft). 

2  y 

P  =  /?sat - COS0  (11) 

r 

Apparent  contact  angle  on  a  rough  porous  surface,  such  as  a 
textile  surface,  is  expressed  using  the  true  contact  angle  ft  [16]. 

cos  0  =  /s  COS  ft  -  fo  (12) 

where  /s  and  fQ  are  fractions  of  solid  and  open  areas  of  the 
hydrophobic  microlayer,  respectively.  If  the  true  contact  angle 
ft,  is  greater  than  90°,  the  apparent  contact  angle  ft  is  increased. 
The  apparent  contact  angle  ft  is  estimated  to  be  between  150° 
and  160°  (under  the  condition  of  a  true  contact  angle  range 
between  120°  and  140°,  and  a  solid  fraction /s  =  0.3).  The  pres¬ 
sure  of  water  in  the  catalyst  layer  is  determined  by  the  maximum 
pore  size  of  the  GDL.  The  increase  of  water  pressure  is  esti¬ 
mated  to  be  approximately  0.08  bar,  using  the  surface  tension 
y  =  0.063  Nm"1  and  a  pore  radius  r=30[xm.  A  pore  radius 


Fig.  4.  Capillary  effect  on  the  net  water  drag  a,  as  a  function  of  cathode  humid¬ 
ified  temperature;  (□)  denotes  the  calculation  without  capillary  effect  and  (■) 
denotes  the  calculation  with  capillary  effect.  (FC  mode,  Cell  =  80  °C,  Ca  =  60  °C, 
£f  =  1.43,  ft  =  2);  (•,  O)  are  the  same  as  shown  in  Fig.  1. 


of  r=30pim  is  applied  because  the  diameters  of  90%  of  the 
pores  in  the  GDL  are  in  the  range  less  than  60  pm.  Net  drags 
were  computed  again,  substituting  Eq.  (1 1)  for  the  partial  pres¬ 
sure  p ,  in  Eq.  (4).  The  results  are  shown  in  Fig.  4.  The  newly 
computed  net  water  drags  (solid  squares)  are  approximately 
0.02-0.07  lower  than  the  corresponding  previous  data  (empty 
squares).  Although  a  discrepancy  still  exists,  the  coincidence 
between  the  experimental  and  computed  data  was  improved. 
This  suggests  a  contribution  of  capillary  effect  (suitability  of 
hydrophobic  effect)  to  the  net  drag. 

The  water  content  in  the  cathode  gas  increases  along  the  flow 
channel  by  the  accumulation  of  generated  water  and  net  water 
drag.  Therefore,  the  water  vapor  in  the  cathode  gas  may  reach 
saturation  pressure  at  a  downstream  position.  The  net  water  drag 
depends  on  two  factors  given  in  Eq.  (2),  i.e.  the  electro-osmotic 
drag  coefficient  n d,  and  the  difference  in  activity  between  the 
cathode  and  the  anode.  Fig.  5  shows  the  net  drag  along  the  flow 
channel  for  various  humidity  of  the  inlet  anode  gas  at  relatively 
low  cathode  humidity,  ft c  =  64°C,  where  the  humidity  ftc,  is 
expressed  using  the  dew  point  of  the  inlet  cathode  gas.  Fractional 
distance  refers  to  the  normalized  distance  from  the  inlet  to  the 
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Fig.  5.  Change  of  net  water  drag  a,  along  a  cathode  flow  channel  for  vari¬ 
ous  humidity  of  the  anode  inlet  (Cell  =  70  °C,  Cc  =62°C).  The  horizontal  axis 
is  the  fractional  distance  along  the  cathode  flow  channel.  (♦)  Ca  =  70°C;  (□) 
Ca  =  65°C;(0%da  =  52°C. 
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Fig.  6.  Net  water  drag  a,  as  a  function  of  cathode  humidified  temperature.  Param¬ 
eters:  (a)  tceii  =  80  °C,  /da  =  80  °C  (•),  60  °C(0);  (b)  Ceil  =  VO  °C,  /da  =  70  °C  (♦), 
52  °C  (0). 

exit.  For  most  of  the  cell  area  from  the  cathode  inlet,  the  net 
drags  do  not  show  any  significant  change  for  a  wide  range  of 
inlet  humidity.  Near  the  cathode  exit  (between  the  normalized 
distance  0.7  and  0.9),  the  net  drag  is  widely  distributed,  from 
0.55  to  —0.34,  for  a  wet  anode  gas  (/da  =  70  °C)  to  a  dry  anode 
gas  (/da  =  52°C).  The  average  net  drag  ranges  between  0.076 
and  0.275,  as  shown  by  the  dashed  line  in  Fig.  5.  This  indicates 
that  the  change  of  net  drag  shown  in  Fig.  1  depends  mainly  on 
the  humidity  near  the  cathode  exit. 

In  the  previous  experiments  (Figs.  1  and  2),  the  net  water  drag 
was  measured  under  a  relatively  low  humidity  anode  gas  con¬ 
dition.  The  net  water  drag  for  various  humidity  (Uc)  of  the  inlet 
cathode  gas  was  calculated,  and  the  results  are  shown  in  Fig.  6 
for  wet  anode  gas  (/da  =  80  °C)  and  dry  anode  gas  (/da  =  60  °C). 
For  a  wet  anode  gas  condition,  the  value  of  net  water  drag  is 
larger  than  that  for  a  dry  anode  gas  condition.  However,  the 
computed  net  drag  for  a  wet  anode  gas  and  dry  cathode  gas  con¬ 
dition  is  approximately  0.6,  which  is  far  less  than  the  measured 
value  of  1.0  at  saturation  [15].  This  difference  is  due  to  back 
diffusion. 

In  order  to  compare  the  present  model  with  the  reported 
experimental  results  [1],  computations  for  different  fuel  cell 
operation  conditions  were  made.  The  conditions  used  for  com¬ 
putation  were  as  follows;  current  density  0.4  A  cm-2,  pressure 
1.5  bar,  cell  temperature  60  °C  and  stoichiometry  of  1.5  and  2 
for  the  anode  and  cathode,  respectively.  The  humidity  of  the 
cathode  is  changed  under  two  different  anode  humidity  (wet 
and  almost  dry)  conditions.  The  computed  results  are  shown 
in  Fig.  7,  with  the  experimental  data,  and  they  are  in  good 
agreement. 


3.2.2.  Sustainability  of  operation  under  the  dry  condition 
When  the  cathode  gas  at  the  channel  is  dry,  the  flux  of  water 
vapor  from  the  cathode  electrode/membrane  interface  to  the 
cathode  gas  channel  may  exceed  the  rate  of  water  production 
and  water  transported  from  the  anode.  This  eventually  causes 
drying  of  the  cathode  catalyst  layer  and  membrane.  However, 
the  performance  of  the  fuel  cell  has  been  shown  to  be  stable  up 
to  an  operating  cell  temperature  of  60  °C,  even  under  a  dry  air 
condition  [17]. 
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Fig.  7.  Comparison  of  experimental  data  with  computed  results  from  the  present 
model.  Experimental  data  are  from  Janssen  [1]  and  is  marked  as  (■).  The  cal¬ 
culation  was  performed  under  the  same  conditions  as  described  in  [1].  Cell 
temperature  and  operating  pressure  are  60  °C  and  1.5  bar,  respectively.  (□) 
/da  =  60  °C;  (O)  /da  =  32  °C.  AD  and  CW  represent  anode/dry  and  cathode/wet 
(/dc  =  60  °C)conditions,  respectively. 


This  result  suggests  the  existence  of  a  mechanism  to  maintain 
minimum  water  content  in  the  membrane  and  the  catalyst  layer 
under  dry  gas  conditions.  A  liquid-vapor  water  interface  with  a 
small  meniscus  exists  in  porous  media.  The  pressure  difference 
across  the  interface  causes  a  change  in  the  free  energy,  which  is 
expressed  by  the  following  Kelvin  equation  [18]: 


P  w 

- =  exp 

Psat 


(13) 


where  pw  is  the  vapor  pressure  above  the  meniscus,  and  pSSLt 
is  the  vapor  pressure  at  the  free  surface.  This  shows  that  the 
equilibrium  vapor  pressure  of  the  meniscus  decreases  from  that 
of  the  free  surface.  Under  a  dry  or  unsaturated  gas  condition, 
the  water  vapor  pressure  at  the  electrode/membrane  interface 
is  distributed  along  the  flow  direction,  and  is  less  than  the  sat¬ 
uration  pressure  for  the  cell  temperature.  Under  the  basic  cell 
operation  conditions  (Uc  =  65  °C),  the  water  vapor  pressure  at 
the  electrode/membrane  interface  is  distributed  from  28.2  to 
31.1  kPa  and  increases  towards  the  exit,  as  shown  in  Fig.  8 
(solid  square).  The  difference  in  vapor  pressure  can  be  explained 
by  the  decrease  in  the  equilibrium  vapor  pressure,  as  shown 
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Fig.  8 .  Change  of  water  vapor  pressure  along  a  cathode  flow  channel  for  different 
humidity  of  the  cathode  inlet  (/ceii  =  70°C,  /da  =  65°C);  (■)  /dc  =  65°C;  (□) 
/dc  =  60°C. 
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in  Eq.  (13).  At  a  fractional  distance  of  0.1,  pores  smaller  than 
45.4  nm  hold  and  keep  water.  For  the  lower  humidity  condition 
(fdc  -  60  °C,  open  square),  pores  smaller  than  14.1  nm  hold  and 
keep  water.  Although  under  dry  and  unsaturated  gas  conditions, 
liquid  water  contained  in  relatively  large  pores  evaporates,  liquid 
water  still  exists  in  micro  pores  and  this  explains  the  experimen¬ 
tal  data  for  stable  cell  performance  observed  even  under  dry 
conditions. 

4.  Conclusion 

Net  drag  was  measured  for  H2  and  FC  mode  operations  under 
various  humidity  conditions.  In  contrast  to  FC  mode  opera¬ 
tion,  flooding  in  the  cathode  GDF  is  not  expected  to  occur 
for  H2  mode  operation.  In  addition  to  these  measurements,  a 
two-dimensional  water  transport  model  for  a  proton-exchange 
membrane  fuel  cell  was  developed.  Net  water  transport  from  the 
anode  to  the  cathode  varies  along  the  flow  channel,  depending  on 
the  net  drag  and  current  density  distribution.  The  distribution  of 
current  density  along  the  flow  channel  was  calculated  using  the 
membrane  ionic  conductivity  and  over-potential  under  a  condi¬ 
tion  of  uniform  cell  voltage. 

The  measured  net  drag  was  compared  with  the  results  com¬ 
puted  using  the  two-dimensional  cell  model.  Although  the  cal¬ 
culated  results  are  slightly  higher  than  the  measured  values,  they 
followed  similar  trends.  However,  there  was  a  larger  discrepancy 
for  the  FC  mode  operation,  and  this  suggests  that  a  resistance 
of  some  kind  exists  for  water  transport  to  the  cathode  channel 
from  the  cathode/electrode  interface.  The  hydrophobic  micro¬ 
layer  that  is  applied  to  the  interface  of  the  catalyst  layer/GDF 
is  considered  to  be  one  candidate  for  the  resistance  to  water 
transport,  and  the  effect  on  the  net  water  drag  was  estimated. 
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